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The process of radiation induced electron capture by protons or deuterons producing new ultra 
low momentum neutrons and neutrinos may be theoretically described within the standard field 
theoretical model of electroweak interactions. For protons or deuterons in the neighborhoods of 
surfaces of condensed matter metallic hydride cathodes, such conversions are determined in part 
by the collective plasma modes of the participating charged particles, e.g. electrons and protons 
or deuterons. The radiation energy required for such low energy nuclear reactions may be supplied 
by the applied voltage required to push a strong charged current across a metallic hydride surface 
employed as a cathode within a chemical cell. The electroweak rates of the resulting ultra low 
momentum neutron production are computed from these considerations. 

PACS numbers: 12.15.Ji, 23.20.Nx, 23.40.Bw, 24.10.Jv, 25.30.-c 



I. INTRODUCTION 

Excess heats of reaction have often been observed to 
be generated in the metallic hydride cathodes of certain 
electrolytic chemical cells. The conditions required for 
such observations include high electronic current densi- 
ties passing through the cathode surface as well as high 
packing fractions of hydrogen or deuterium atoms within 
the metal. Also directly observed in such chemical cells 
arc nuclear transmutations into elements not originally 
present prior to running a current through and/or prior 
to applying a LASER light beam incident to the cath- 
ode surface [J H, [H, 0, H, @- It seems unlikely that the 
direct cold fusion of two deuterons can be a requirement 
to explain at least many of such observations [7| because 
in many of these experiments, deuterons were initially 
absent. For simplicity of presentation, we consider "light 
water" chemical cells in which deuterons are not to any 
appreciable degree present before the occurrence of heat 
producing nuclear transmutations. 

Nuclear transmutations in the work which follows are 
attributed to the creation and absorption of ultra low 
momentum neutrons as well as related production of 
neutrinos. Although other workers jioi hjj have previ- 
ously speculated on a central role for neutrons in such 
transmutations, they were unable to articulate a physi- 
cally plausible mechanism that could explain high rates 
of neutron production under the stated experimental con- 
ditions. By contrast, in this work electrons are captured 
by protons all located in collectively oscillating "patches" 
on the metallic surface. Electrons are captured by pro- 
tons all located in collectively oscillating "patches" on 
the metallic surface. Since the energy threshold for such 
a reaction is 

Qin ~ {M n - (M p + to)} c 2 « 0.78233 MeV, (1) 
one requires a significant amount of initial collective radi- 



ation energy to induce the proton into neutron conversion 



(radiation energy) 



P 



(2) 



The radiation energy may be present at least in part due 
the power absorbed at the surface of the cathode. If V de- 
notes the voltage difference between the metallic hydride 
and the electrolyte and if J denotes the electrical current 
per unit area into the cathode from the electrolyte, then 
the power per unit cathode surface area dissipated into 
infrared heat radiation is evidently 



V = VJ = eV$, 



(3) 



wherein $ is the flux per unit area of electrons exiting 
the cathode into the electrolyte. Typical metallic hy- 
dride cathodes will exhibit soft surface photon radiation 
in much the same physical manner as a "hot wire" in a 
light bulb radiates light. For the case of chemical cell 
cathodes, there will be a frequency upward conversion 
from virtually DC cathode currents and voltages up to 
infrared frequency radiation. Such an upward frequency 
conversion requires high order electromagnetic interac- 
tions between electrons, protons and photons. 

The purpose of this work is to estimate the total rates 
of the reaction Eq.@ in certain metallic hydride cath- 
odes. The lowest order vacuum Feynman diagram for 
the proton to neutron conversion is shown in FIG.[TJ For 
the case of the reactions in metallic hydrides, one must 
include radiative corrections to FIG. [1] to very high order 
in the quantum electrodynamic coupling strength 
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7.2973526 x 10 



-3 



(4) 



The W^-coupling in terms of the weak rotation angle Ow 
will be taken to lowest order in 



aw 
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Charge conversion reactions are weak due to the large 
mass My/ of the W^. The Fermi interaction con- 
stant, scaled by either the proton or electron masses, is 
determined by 



G, 



he 



G F Ml 

- v « 1.02682 x 10~ 5 , 

he 

G F m 2 



he 



3.04563 x 10~ 12 . 



(6) 



In the work which follows, it will be shown how the 
weak proton to ultra low momentum neutron conversions 
on metallic hydride surfaces may proceed at appreciable 
rates in spite of the small size of the Fermi weak coupling 
strength. 

An order of magnitude estimate can already be derived 
from a four fermion weak interaction model presuming 
a previously discussed [HJ electron mass renormalization 
m — ► m = (3m due to strong local radiation fields. Sur- 
face electromagnetic modes excited by large cathode cur- 
rents can add energy to a bare electron state e~ yielding 
a mass renormalized heavy electron state e~ , with 



m = (3m. 



(J) 



The threshold value for the renormalized electron mass 
which allows for the reaction Eqs.([T]) and ^ is 



(3 > (3 ~ 2.531. 



(8) 



For a given heavy electron- proton pair (e _ p + ), the tran- 
sition rate into a neutron-neutrino pair may be estimated 
in the Fermi theory by 



G Fm 2\ 2 /_2 



he 
9 x 10 



-j- ) (P-PoY, 



T {s - p+) ^ n+Ve ~ 7 x 10- 4 Hz x (/3 - A,) 2 , (9) 

If there are n-i ~ 10 16 /cm 2 such (e~p + ) pairs per unit 
surface area within several atomic layers below the cath- 
ode surface, then the neutron production rate per unit 
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FIG. 1: A low order diagram for e~ + p — > n + v e in the 
vacuum is exhibited. In condensed matter metallic hydrides, 
the amplitude must include radiative corrections to very high 
order in a. 



surface area per unit time may be estimated by 



n7 2 i=i ri2T^ p +y 



>n+z/ e > 



TU 2 



10 13 Hz 



cur 



x {(3-(3 ) 2 . 



(10) 



Significantly above threshold, say (3 ~ 2(3$ ~ 5, the esti- 
mated rate ~ 10 13 Hz/cm 2 is sufficiently large so as to 
explain observed nuclear transmutations in chemical cells 
in terms of weak interaction transitions of (e~p + ) pairs 
into pairs into neutrons and neutrinos and the subse- 
quent absorption of these ultra low momentum neutrons 
by local nuclei. 

It is worthy of note that the total cross section for the 
scattering of neutrons with momentum p may be writ- 
ten a tot = {^h/p)^mA{p) wherein the forward scatter- 
ing amplitude for neutrons is »4(p)[l3[. In the ultra- low 
momentum limit p — ► 0, the cross section for neutron 
absorption associated with a complex scattering length 
b = limp_^o 3m^4(p) formally diverges. For a finite but 
large neutron wavelength A, the mean free neutron ab- 
sorption path length A corresponding to n a absorbers per 
unit volume obeys 



A" 



n a <Jtot ~ 2nXb. 



(11) 



Numerically, for n a ~ 10 22 /cm 3 neutron absorbers per 
unit volume with an imaginary part of the scattering 
length b ~ 10~ 13 cm and with ultra-low momentum 
neutrons formed with a wavelength of A ~ 10~ 3 cm, a 
neutron will move on a length scale of A 



10~ 3 
- 10" 



before being absorbed. The externally detectible neu- 
tron flux into the laboratory from the cathode is thereby 
negligible [T^ • 

Similarly, there is a strong suppression of gamma-ray 
emission due to the absorption of such rays by the heavy 
surface electrons [l3j]. The mean free path of a photon 
in a metallic condensed matter system is related to the 
conductivity er; employing Maxwell's equations and the 
equations 



— = Rvac<J ~ 4a 



(!) 



1/3 



,2/3/ 



(12) 



wherein n is the density per unit volume of heavy elec- 
trons on the cathode surface and I is the heavy elec- 
tron mean free path. For hard photon energies in the 
range 0.5 Megavolt < (hco 7 /e) < 10 Megavolt we esti- 
mate n 2 / 3 ~ 10 15 /cm 2 and I ~ 10 -6 cm leading to the 
short gamma-ray mean free path L 1 ~ 3.4 x 10~ 8 cm. In 
estimating such a small mean free path for hard gamma 
photon, let us remind the reader of the inadequate na- 
ture of single electron single photon scattering when mak- 
ing estimates of photon mean free paths. For example, 
a single optical photon scatters off a single single non- 
relativistic electron with a Thompson cross section. For 
many electrons, if the Thompson cross section were ap- 
plied then one might conclude that metals were transpar- 
ent in the optical regime. The proper estimates of photon 
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FIG. 2: The four fermion vertex for e~ + p — > n + v e in the 
vacuum is exhibited. In the large Mw limit, the Feynman 
diagram of FIG. Q] collapses into the above Feynman diagram. 
In condensed matter metallic hydrides, the resulting effective 
W ± fields are defined in Eqs.dEJ and (J39J) . 

absorption rates in condensed matter employs the elec- 
trical conductivity. 

In SecHIl an exact expression is derived for the emis- 
sion rate w per unit time per unit volume for creating 
neutrinos. It is then argued, purely on the basis of con- 
servation laws, that vj also represents the rate per unit 
time per unit volume of neutron production. The rate 
w, in Sec lIIH is expressed in terms of composite fields 
consisting of charged electrons and opposite charged W- 
bosons. The effective IF-bosons for condensed matter 
systems may be written to a sufficient degree of accuracy 
in terms of Fermi weak interaction currents 

^ = C {ipnlnigv - 9AlB)lpp) i 

X » = c ($ P Jn(gv - £U7s)V>n) , (13) 

wherein the Dirac matrices are defined in Sec|IIJ ip p and 
ij) n represent, respectively, the proton and neutron Dirac 
fields and the vector and axial vector coupling strengths 
are determined by 

A= — re 1.2695, 

9V 

cos 6 C = g v « 0.9742, (14) 

wherein 9c is a strong interaction quark rotation angle. 
In Sec lIVl the electron fields as renormalized by metallic 
hydride surface radiation are explored and the effective 
mass renormalization in Eq.((7]) is established. In SecfVl 
we consider the coupled electron and proton oscillations 
near the surface of a metallic hydride. 

In Sec lVII the nature of the neutron production is dis- 
cussed in terms of isotopic spin waves. In the limit 
in which the protons and neutrons are non-relativistic, 
one may view the proton and neutron as different iso- 
topic spin states of a nucleon[3] with the charged pro- 
ton having an isotopic spin +1/2 and with the neutron 
having an isotopic spin —1/2. If n(r) and p(r) repre- 
sent, respectively, the two (real) spin component fields 
for non-relativistic neutrons and protons, then the oper- 
ator isotopic spin density T(r) = (7i(r), 72 (r),^^)} of 



the many body neutron-proton states may be written 



T-i (r ) 


~ 2 ^ 


lb yi )p\L )) , 


Ta(r) 


= 2 ( nt ( r M r ) - 


■p\r)n{r)) , 


Ts(r) 


= \ (PH r )p( r ) - 


n f (r)n(r)) , 


T±(r) 


= Ti(r)±*T 2 (r) 





(15) 



In the non-relativistic limit, these isotopic spin operators 
determine the time-component of Fermi weak interaction 
currents in Eq. (|13p via 

ZT0 (r)wc5yT ± (r) . (16) 

The remainder of the non-relativistic weak currents are 
of the Gamow- Teller variety[l5| and require the true spin 
as well as isotopic spin version of Eq. (|16p : i.e. with S = 
cr/2 as the Fermion spin matrices, the combined spin and 
isotopic spin operator densities are 

S ltj (v) = ( P Hr)S 3 n(r)+nHr)S 3 p(r)), 

5 2j (r) = i(n\r)S jP (r)-pH*)SMr)), 

5 3j (r) = (pt( r )5 i: p(r)-nt(r)^n(r)), 

^^(r) = S 1<j (r)±iS 2 , j (r). (17) 

In the non-relativistic limit for the protons and neutrons, 
the spatial components of the weak interaction currents 
in Eq.((T5)| are 

I^(r)«-c^(r). (18) 

Altogether, in the nucleon non-relativistic limit 

I T p w c (-g A Sf, ~g A St, -gASf.gv^) . (19) 

The isotopic formalism describes the neutron creation as 
a surface isotopic spin wave. Out of many oscillating 
protons in a surface patch, only one of these protons will 
convert into a neutron. However, one must superimpose 
charge conversion amplitudes over all of the possibly con- 
verted protons in the patch. This describes an isotopic 
spin wave localized in the patch with wavelength 
The wavelength in turn describes the ultra low momen- 
tum p ~ hk of the produced neutron. Finally in the 
concluding Sec lVIIl further numerical estimates will be 
made concerning the weak interaction production rate of 
such neutrons. 



II. NEUTRINO SOURCES 

The conventions here employed are as follows: The 
Lorentz metric r rf ,v has the signature (+, +, +, — ) so that 
the Dirac matrix algebra may be written 

7 ^ 7 " = - r f» - ia^ wherein = (20) 
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The chiral matrix 75 is defined with the antisymmetric 
symbol signature e 12 3o = +1 employing 



1 
4! 



11 v X a 



and chiral projection matrices are thereby 



1 



^± = 2(^75). 



(21) 



(22) 



Further algebraic matrix identities of use in the work 
below, such as 

h \^v = ie X„<TU _ + ^Xa^u _ r ^r T ^ } ( 2 3) 



all follow from Eqs.(|20]), d2Tj) and j22|). 

The average flux of left handed electron neutrinos (pre- 
sumed massless) is determined by 



T»{x) = c(v(x)^P+v(x)) . 



(24) 



Initial state averaging in Eg. (f2~4")l is with respect to a 
chemical cell density matrix 



>...)=Trp(...), 



(25) 



The mean number of neutrinos created per unit time per 
unit volume may be computed from the four divergence 
of the neutrino flux; i.e. 



w{x) = d^(x). 



(26) 



Let us now argue, purely from standard model conserva- 
tion laws, that m is also the mean number of neutrons 
created per unit time per unit volume within the metallic 
hydride cathode in a chemical cell. 

If a neutrino is created, then lepton number conserva- 
tion dictates that an electron had to be destroyed. If an 
electron is destroyed, then charge conservation dictates 
that a proton had to be destroyed. If a proton is de- 
stroyed, then baryon number conservation dictates that 
a neutron had to be created. Thus, the rate of neutrino 
creation must be equal to the rate of neutron creation. It 
is theoretically simpler to keep track of neutrino creation 
within the cathode. 

The neutrino sinks and sources, respectively fj and 77, 
are defined by that part of the standard model action 
which destroy and create neutrinos; i.e. 



Siii 



(r](x)v(x) + v{x)rj(x)) d 4 x. 



(27) 



The neutrino field equations are thereby 

- i"f^d^v(x) = r](x), 
id tll v(x)^ = fj(x). 



Eqs. (|24|) . ((26)) and (|28[) imply the neutrino creation rate 
per unit time per unit volume at space-time point x has 
the form 



j(x) — 2c3m {rj(x)P + v(x)) 



(29) 



Introducing the retarded massless Dirac Green's func- 
tion, 

-i^d fl S(x-y)=6(x-y), (30) 
allows us to solve the neutrino field Eqs. (|2"5)) in the form 

wherein Vi n {x) represents the asymptotic incoming neu- 
trino field. The assumption of zero initial background 
neutrinos is equivalent to the mathematical statement 
that the neutrino destruction operator v1[ n {x) \T) = for 
the initial states in Eq. ([231) . In such a case, Eqs. (|2T)|) and 
(EH) imply 

w{x) = 2c3m / (fj(x)P+S(x - y)r,(y)) d 4 y. (32) 



The retarded massless Dirac Green's function may be 
found by looking for a solution of Eq. (13"01) of the form 

S{x-y)=i'fd ll A(x-y). (33) 

From Eas.([5D| and Q it follows that 

- d^A(x -y) = S(x - y). (34) 



The retarded solution to Eq. ([34|) requires the step func- 
tion 



In detail 



tf(x - y) = l if x° > y°, 
®(x -y) = if x° < y°- 



A{x-y) = § Xr jl 5({x-y?) 

Z7T 



(35) 



(36) 



Eqs.(32]) and (02]) imply 

w(x) = 2c5Re / {fj{x)P + ^r,{y)) d^A(x - y)d A Vl 
w(x) = 2c5Re / A(x - y) {fj(x)P + ^r,(y)) d 4 y. (37) 



(28) 



The neutron production rate w per unit time per unit 
volume can thus be computed in terms of the neutrino 
sinks fj and sources r\. 



III. COMPOSITE CHARGED FIELDS 

The neutrino sinks and sources of interest in this work 
can be written in terms of the composite fields of charged 
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electrons and charged effective condensed matter W - 
bosons; i.e. 

r)(y) = ±gfW+{y)P+1>{y) , 
fj(x) = ±=${x)P^ x W^{x) , (38) 
in which ip and ip are the Dirac electron fields and 
WM = (^-) #(y), 

Wx(x)=(^^)l^{x) . (39) 

The weak interaction proton-neutron charged conver- 
sion currents are defined in Eq. (|13[) . In Feynman 
diagram language, the amplitude pictured in FIG[T] has 
been replaced via a field current identity of the Fermi 
four field point interaction in FIGH Eas.([3"7|) and l[3"5j) 
imply 

w{x) = c 5Re J g{x, y)A(x - y)a A y, (40) 

wherein 

G{x,y) = 

h X^U ( W -( x )$( x ) 7ud ^p + ^y) W +( y ))) . (41) 

The neutron production rate w per unit time per unit 
volume implicit in Eqs. (f3"6")) . (|4"0")) and (|4"Tj) may be con- 
sidered to be exact. 

IV. ELECTRON MASS RENORMALIZATION 

When the reacting charged particles e~ + p + — > n + v e 
are in the presence of surface plasmon radiation, then the 
external charged lines (incoming wave functions) must 
include the radiation fields to a hi gh o rder in the quantum 
electrodynamic coupling strength 16] a. The situation is 
shown in FIGJ3] To see what is involved, recall how one 
calculates the density of states for two particles incoming 
and two particles outgoing: 

Case I: The Vacuum The density of states may be written 
as the four momentum conservation law, 

J e i{s 0+ +s --s 0n -s „)/h d i x _ 

{2irh) 4 S( P+ +p_- Pn -p„), (42) 

wherein p + , p_, p n and p v represent, respectively, the 
four momenta of the proton, electron, neutron and neu- 
trino. For a particle of mass m in the vacuum, the action 
Sq(x) obeys the Hamilton- Jacobi equation 

d fl S (x)d fl S (x) +m 2 c 2 = 0, 

S (x)=p-x= Pfl x>*. (43) 




Surface Radiation 



FIG. 3: In the presence of electromagnetic surface radiation, 
the charged particles in weak interaction must be described 
by wave functions which include to high order in a the effects 
of the electromagnetic fields. For the reaction at hand, both 
the proton and the electron react to surface radiation. The 
resulting mass renormalization is stronger for the electronic 
degrees of freedom than for the proton degrees of freedom. 
The density of states including radiation is computed em- 
ploying Eqs.((45} and (146) . 

Case II: Radiation If the reaction takes place in the pres- 
ence of electromagnetic radiation, 

= - d v A^ (44) 

then the density of states conservation of four momenta 
must also include the electromagnetic radiation contri- 
bution; i.e. 

Me [ e i(s + +s—s n -s„)/K d 4 x = 

(2TTh) i S{ P+ ,p-,p n ,p I/ ), (45) 

wherein, for a charged particle, the Hamilton- Jacobi 
equation reads [l7j 

mvn(x) = d M S(x) - -Apfa), 

v^(x)v^(x) +c 2 = 0. (46) 

Therefore, in the density of states Eq. (|4"5|) including ra- 
diation the full solution of the Hamilton-Jacobi equation 
must be solved for all of the charged particles in the inter- 
action. This constitutes the physical difference between 
the diagrams in the vacuum shown in FIG (Hand includ- 
ing radiation shown in FIG [3] Under a gauge transfor- 
mation A^ — > A^ + d^X-i the Hamilton-Jacobi Eq. (|4"6"|) 
for a charged particle implies a change in the action 
S — ► S + ex/c. The renormalized charged particle wave 
function thereby exhibits the expected gauge transfor- 
mation rule ip ~~ * ip exp(ex/^c) making the complete am- 
plitude from FIG [3] gauge invariant. The electron gauge 
contribution to the phase is canceled by the proton gauge 
contribution to the phase since the two particles are op- 
positely charged. 
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The mass renormalization may be understood by av- 
eraging the local momentum p^ = d^S over local space 
time regions. Presuming p^A^ = we have on average 
that 

-p~pT^in 2 c 2 + {^j A^A^ = m 2 c 2 . (47) 

The mass renormalization parameter in Eq.© is then 
given by 







1 



(48) 



Since the electron mass is much less than the proton 
mass, m <C M p , the main effects on low energy nuclear 
reactions are due to the mass renormalization of the sur- 
face electrons [T^|. 

From the viewpoint of classical physics, the Lorentz 
force on a charge equation of motion, 

mc— — = eF^ v„, 
or 



(49) 



\y/l- (V/C) 2 ' y/l-(v/c)*J ' 

is reduced to first order via the Hamilton- Jacobi Eq.(|46jl. 
according to 



~dr 



(50) 



From the viewpoint of quantum mechanics, there is a one 
to one correspondence between quantum solutions of the 
Dirac equation and the classical solutions of Hamilton- 
Jacobi equation. In detail, the Dirac equation in an ex- 
ternal radiation field, 



(x) I ip(x) + map(x) = 0, (51) 



may be subject to a non-linear gauge transformation 
employing the solution to the classical Hamilton- Jacobi 
equation, 

?P( X ) = e lS{x)/h ^{x). (52) 

The resulting radiation renormalized wave function obeys 

7 M (-ihd^ + mvn(x)) *(x) + mcf(i) = 0. (53) 

It is worthy of note in the quasi-classical limit H — > 
that the solution to the charged particle wave Eq. (f53"f is 
reduced to algebra. 

V. PROTON OSCILLATIONS 

Thermal neutron scattering from hydrogen [TH [l||, 
metallic hydrides [2(| EE H^i and even from protons 
embedded in proteins [24j has been of considerable recent 



theoretical and experimental interest [25j,|26|. Employing 
the total cross section per unit volume as the extinction 
coefficient, 



hfnt — 



Ctot 



(54) 



the intensity of a neutron beam when passed through a 
sample of thickness d obeys 



1(d) = J(0)expHwO- 



(55) 



The neutron-proton cross section is quite large since there 
is a very shallow bound state (the deuteron). For high 
densities of protons, the neutron-proton scattering will 
normally dominate the condensed matter scattering of 
thermal neutrons. Let the momentum transfer and en- 
ergy transfer to the neutron be defined, respectively, as 



P/ 



and hw 



(56) 



Decomposing the neutron scattering into various momc- 
menta and energies yields 



Af ( Pf\ ddi^ 



dflfdef HV \p. 



dfl 



f 



(57) 



wherein da^f is the differential cross section for a single 
neutron to scatter off a single proton and the dynamical 
proton form factor for Af protons is defined as 



AfS(Q,t 



E 



OO 1 1 



3=1 

Fj(Q,t) = / e -iQ-RiC*) e iQ-R 3 (o) 



2tt 



(58) 



In what follows we shall consider the self diffusion for a 
fixed proton at position Rj(t) = R(i) and write 



e -iQ-R(t) e iQ-R(0) 



S(Q,u)e~ ibri: du. (59) 



Under the assumption that the proton motions are non- 
relativistic, one easily shows for the many body Hamito- 
nian TL that 



Hq = e- iQR He iQR = H + hV-Q 
wherein the velocity of the proton is 



h 2 Q 2 

2M 



Thus 



F(Q,t) 
F(Q,t) 
F(Q,t) 



V = -i — V. 

M 



e -iQ-R(t) e »Q-R(0) 



JHt -iQ-R — iHt _iQ-R 



(60) 



(61) 



■ R ) 



_ / e mt/h e -in Q t/h 



(62) 
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wherein "+" indicates time ordering, 



B. Recoilless Fraction 



HQ 2 
2M 



and V(s) = e ins/h Ye -iHs/h^ ^ 



From Eqs. ([55|) . (|6"2")l and (|63[) one derives the sum rules 



S(Q,ui)dcj = 1, 

/oo 
-OO 



(64) 



which represent, respectively, probability normalization 
and the fact that the mean recoil energy is the same as 
would have been computed for a single free proton. We 
can now discuss the experimental measurements which 
make proton oscillations unusual. 



A. Deep Inelastic Neutron Scattering 

For energetic neutrons with relatively high momentum 
transfer, it was estimated that the neutron scattering 
event lasted on a time scale of an attosecond. It was 
thought that the velocity could not appreciably change 
in so short a time. Therefore one might try the impulse 
approximation wherein Eq. (f6"2")) could be estimated by 
the t — > limit 



^impui S e(Q,t) = e^ Q *<e^ Q - vt ) 



S lmpulse (Q, w) = J S(u - Q • V - uj Q )p(V)d 3 V, (65) 



wherein p(V)d 3 V is the probability that the proton ve- 
locity is in the range V £ d 3 V. Eq. ([65|) is at the heart of 
the theoretical analysis of so-called "deep inelastic" neu- 
tron scattering. Note that the impulse approximation 
obeys the sum rules in Eq. (|64p . It came as somewhat of 
a shock that for protons inside liquid hydrogen, metallic 
hydrides and proteins, the impulse approximation does 
not work. At issue was the "violation" of the probability 
sum rule 



L 



S{Q,u)duj < 1. 



experimental 



(66) 



Since Eqs. ([54"|) and (|55|) were experimentally verified, the 
only reason for the "experimental probability loss" must 
have been that the integral over experimental high fre- 
quency data was missing a physical low frequency regime 
which persists and contributes to the probability integral. 
Such low frequency modes must be collective. 



Let us at first assume the cluster decomposition prop- 
erty of F(Q,t); i.e. 

lim F(Q,t) = Foa(Q), 



(67) 



t Hm(e- iQ - R We iQ ' R (°)) = |(e 4Q ' R )|- 
F 00 (Q) = |(e 4Q - R )| 



If cluster decomposition holds true, then F OQ (Q) repre- 
sents the recoilless fraction in the sense that 



5(Q, u) = F 00 (Q)S(uj) + S(Q, to) 

S(Q,a;)dw = l-F 0O (Q)<l. (68) 



Since the second of Eqs. (f6"8")) yield the fraction that re- 
coils, the experimental puzzle in Eq. (|66[) is then resolved. 

The recoilless fraction ^^(Q) for the proton is the 
same factor as the recoilless fraction F oa (Q) that appears 
in the Mossbauer effect for the gamma ray decay in a 
heavy nucleus. In the case of the Mossbauer effect, the 
recoil of the gamma emission is taken up by the crystal 
as a whole. For the case of deep inelastic neutron scatter- 
ing, the proton recoil must be taken up by neighboring 
electrons and other protons in coherent oscillation. Let 
P(u)d 3 u be the probability that the displacement of a 
proton from its equilibrium position u = R — ^R^ is in 
the range u S d 3 u. The recoilless fraction in Eq. (|6"7|) may 
be written as 



FooiQ) 



J P(u)e lQ u d 3 u 



(69) 



The mean square fluctuations in displacement as defined 
by 



uu = J F(u)uurf 3 u, 



(70) 



can be measured via the low momentum transfer limit 

mFoo(Q) 



lim ; ■ = — lim 



Q^o 



(71) 



For higher momentum transfers, the fractal dimensional 
spectator model 24] gives a fairly good representation of 
the recoilless fraction; i.e. 



Foo(Q;£>) - 



IQ-u| : 



D 



(72) 



The harmonic oscillation result is the Gaussian probabil- 
ity Debye- Waller factor 



lim Foo(Q; D) = exp -|Q ■ u 



(73) 
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At elevated temperatures, a small value of D describes 
the fractal dimension of the probability distribution 
P(u)d 3 u. Finally, for the case of collective proton oscilla- 
tions, the delta function in Eq. (|68[) is somewhat idealized, 
In experimental practice, the delta function is broadened 
into an infrared peak consistant in thermal equilibrium 
with the detailed balance condition 



S(— Q, —oj) = e 



_ —Hw/ki 



T S(Q, 



(74) 



C. Proton Displacements 



Neutron scattering experiments [20] on palladium hy- 
dride at moderate momentum transfer clearly indicate a 
sharply defined collective oscillation peak at (Ml/e) w 
60 millivolt. Such a collective proton oscillation at an 
infrared frequency will resonate with electronic surface 
plasmon oscillations of the electrons leading to the local 
breakdown of the Born-Oppcnhcimcr approximation and 
large collective proton oscillation amplitudes. Some the- 
oretical insights into the nature of such oscillations can 
be obtained from sum rules. The mobility tensor of a 
single proton at complex frequency may be defined via 
the Kubo formula 



m{0 = hJo <[ V W' R (°)]> e< ^' 9m(C)>0. (75) 



The mobility tensor obeys the dispersion relation 
-2i(\ f°° 5Rem(w + i0+)duj 



m(C) 



(76) 



As the complex frequency |£| — > oo in the upper half 
plane, the mobility tensor obeys the asymptotic form 



m(C) 



5ftem(w + i0 + )dui 



2i \ f°° 

— r / w 2 5Rem(w + i0+)duj 
"V / Jo 



(77) 



The same asymptotic expansion can be obtained by con- 
tinually integrating Eq. (|75j) by parts yielding the sum 
rules 



3?em(w + i0+)dcj = - ([V, R]) , 
Lu 2 ftem(uj + i0+)duj = i ([A,V]) , (78) 



wherein the proton velocity is V = R and the proton 
acceleration is A = V. Employing the non-relativistic 
proton Eqs.(|6T|). d78|) and Newton's law MA = |e|E, we 
find that 

2 f°° 1 
-/ 3?em(w + i0 + )du = — . 
7T J M 

-J uj 2 ftem{uj + iO+)duj = --y^(VE}. (79) 



Employing Gauss'law at the proton position (diwE) = 
— \e\fi wherein h is the electron density right on top of 
the proton, 



= Vi(R-r t ) =(ft(R W R)) 



fc 



(80) 



the sum rule estimate for the infrared frequency peak in 
the neutron scattering is thereby 



9. z 



/ °° LJ 2 ^e tr m(uj + i0+)doj 
/ °°5fe tr m(u + i0+)doj 
e 2 h 47r ( h 2 
3M _ ~: 



Mm 



(81) 



In Eg. (|8T|) , a is the Bohr radius. In order to understand 
more clearly our estimate of the proton oscillation fre- 
quency O, suppose that the proton was embedded in a 
sphere with charge density — \e\n. If the proton is pulled 
out by a displacement u, then an electric field would exist 



HE 



-MQ 2 u 



(82) 



pushing the proton back to the sphere center. The equa- 
tion of motion Mil — — |e|E = — Mfl 2 u yields the re- 
quired oscillation. One may also note the identity 



IE I 2 = M 2 n 4 '\u 



(83) 



relating the meansquare electric field at the proton posi- 
tion to the mean field displacement fluctuation. Sup- 
pose a mo mentum fo r a relativistic electron given by 
p = mv / -Jl — |v/c| 2 . The energy K of such an elec- 
tron may then be written 



K 2 = m 2 c 4 + |p| 2 c 2 



equivalent to 



K 2 = mV + 



|E| 2 c 2 



^2 4 

m c 



£ 2 



(84) 



(85) 



wherein the electric field scale has been introduced 



mc / Ml 

1 = T [ ~ 



(86) 



Experimentally [20(, for protons absorbed in Pd, one 
finds the soft Boson mode with (hil/e) «6x 10 -2 volts 
wherein 



£ P ePd ~ 1-55 x 10 volts/meter 



(87) 



In this regard, one notes that typical atomic physics elec- 
tric fields in atoms are of that order when located a Bohr 
radius a from an isolated proton; i.e. with the Bohr ra- 
dius 



em amc 



0.5291772108 x 10~ 8 cm, (88) 
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20 



15 - 



CO- 10 




FIG. 4: The predicted electron mass enhacemnt /3 = m/ra 
when protons are absorbed into palladium is plotted as a 
funct ion o f the root mean square proton displacement u where 
u 2 = Id 5 and the Bohr radius a « 0.5291772108 x 10~ 8 cm. 



atomic electric fields are of the order 

M 



£a = 



Ana 2 



5.142206318 x 10 11 volts/meter. (89) 



Employing Eqs.([5T|l. l[55]l and along with the esti- 
mate h w ^(O)) 2 = l/(7ra 3 ) yields 



|E| 2 16 / |u| : 



(90) 



Thus, Eq.([l8"|) may be written^ as 



K 2 



\ 



e 2 |E| 



H^W I (91) 



For example, for protons in palladium Eas.([87|) and ([891 
imply 



which is shown in FIG|4j 



\ 



1 + 19.6 



(92) 



D. Driven Oscillations 

When an electrode has a current density passing into 
the cathode surface, then there is a power per unit area V 



suppled to the surface as in Eq.Q. Under such circum- 
stance the surface temperature rises to a high value. In 
detail, in a non-equilibrium steady state situation, one 
must introduce a "noise temperature" Tq w which de- 
pends on the frequency and wave number of the excita- 
tions which form due to an Ohmic heating rate per unit 
time per unit area V = VJ. The noise temperature en- 
ters into the dynamic form factor via 



S(-Q,-w) = e - huJ / kBT ^ 



S(Q, 



(93) 



which is the non-equilibrium noise temperature version 
of the thermal equilibrium detailed balance Eq. ([74]) . 

Experimentally, the noise temperature is high in the 
infrared w-regime in a small surface domain size d re- 
lated to the wave number via d ~ Q^ 1 . Note that the 
overall temperature over the whole electromagnetic band- 
width need not be very high, but in the infrared regime 
at lj ~ n there should be a sharp peak in the noise tem- 
perature. What has been observed[I3] for deuterons ab- 
sorbed in palladium, are surface hot spots flashing in 
different domains analogous to the flashing of fire-flies 
during a dark evening in the country. The flashing turns 
on and off in apparently random spatial and temporal 
patterns. The hot spots occur when and where regions 
of the rough cathode surface form electromagnetic cavi- 
ties of size and shape as to naturedly support the infrared 
radiation. These occur randomly on rough cathode sur- 
faces. 

During these flash events the nuclear transmutations 
can take place. The hot spots at any given time take 
up only a small fraction of the cathode surface. In- 
creasing the relative fraction of nuclear-active cathode 
surface area which comprises hot spots would propor- 
tionately increase the measured efficiency of neutron 
production on the cathode. The higher the average 
density of hot spots on a given cathode surface over 
time, the greater the amount of excess heat that would 
likely be observed at the cathode device level using 
gross thermal measurement techniques such as calorime- 
try. From the nature of the surface damage on cath- 
odes due to micron-scale hot spots, it is evident that 
the metal actually melts locally during a hot flash pro- 
ducing nuclear transmutations. Distinctive areas of ob- 
vious melting and explosive blow-out cratering features 
that are consistent with the presence of such hot spots 
are very commonly seen in post-experiment SEM im- 
ages of the surfaces of various metallic cathode mate- 
rials that include gold (melting point 1554.9 °C; boil- 
ing point 2963 °C); palladium (melting point 1554.9 °C; 
boiling point 2963 °C)[28|, |29(; titanium (melting point 
1064.18 °C; boiling point 2856 °C)[30], and tungsten 
(melting point 3422 °C; boiling point 5555 °C)^l|. Such 
evidence is consistent with the noise temperature at a 
typical value of the peak Tq w ~5x 10 3 °K correspond- 
ing to a proton displacement of (u/a) ~ 10 within a do- 
main. In accordance with FIG|H the mass enhancement 
is more than sufficent for the neutron production Eqs.|T|) 
and ©. 
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VI. NEUTRONS AND ISOTOPIC SPIN WAVES For steady state production rates, Eq.flM]) reads 



The sources of the neutrinos are inhomogeneous in spa- 
tial regions near the surfaces of cathodes. Also, the neu- 
trinos are so weakly interacting that after emission they 
are virtually unaware of the condensed matter. The neu- 
trino on energy shell phase space Q v = (Q, |Q|) has the 
Lorentz invariant phase space 



dL Q = 



d 3 Q 



2(2tt)3|Q| 



(94) 



Writing the neutrino emission part of Eqs. ([40|) and ([4Tj) 
as the phase space integral 



w{x) 



e iQ-{x-y) h \novQ^ x 



(W~{x)^{x) lv P + 4>{y)W+{y))d A ydL Q 



-Ac 



\ c 4 



(l^{x)^(x)P. 1 ^{y)l+{y))d i ydL Q . (95) 



Under the assumption that the initial proton spins are 
uncorrelated and that the free neutron density is dilute, 
considerable simplification can be made in estimating the 
rather daunting but rigorous Eq. ([95"|) . The estimate for 
the inhomogeneous ultra low momentum neutron pro- 
duction rate per unit volume is 



U 3 ) \ ^ 



(9v + 3si) x 



e «3-(-») (T+(x)$(x)il>(y)T-(y)) d 4 ydL Q ,(96) 



"Re 



wherein Eq. (|19p has been invoked. 

If the neutrons are dilute, then it is sufficient to con- 
sider the creation of a single neutron from a proton, i.e. 
the propagation of the within condensed matter. 
What is left of the heavy boson is merely an iso- 
topic spin wave. There is a superposition of amplitudes 
summed over all the possible protons within a patch 
which may be converted into a neutron. The isotopic 
spin wave creation and annihilation operators in the sur- 
face patch obey 



T^z) «T ± (x)e 



p(cAM)i°/S 



(97) 



with the neutron-proton mass difference determining the 
threshold value of the electron mass m renormalization 
parameter /3 i.e. 



AM = M n — M = p m. 



(98) 



For the creation of a single ultra low momentum neutron 
from non-relativistic protons 



(T + (x)$(x)1>(v)T-(v))=> 
<5(x - y}e -^M)( x «-y°)/n x 

(p^(x)ip(x)ip(y)p(y)) . 



(99) 



hm{r) ^ ^ ( hGp 



(9v + 3.9a) x 



e i(c 2 AM+hc\Q\)t/h x 



(pt( r )^(r)^(r,t)p(r,t)) (cdt)dL Q . (100) 

Explicitly exhibiting the neutrino energy being radiated 
away in Eq.([2]), yields 



hw(r) 1 fhG F 



trie 



2 27r 2 c 2 V 



(9v + 35a) x 



oo />oo 



3?e / / e" 

/O J-oo 



a i(c 2 AM+hc\Q\)t/h x 

(p^(r)ilj(r)i/}(r,t)p(r,t)) (cdt)(u)dw). (101) 



The remaining correlation function in Eq. (|100|) de- 
scribes how an electron which finds itself directly on top 
of a proton propagates in time. The integral over time 
may be written 

/oo 
(pHr)$(r)iP(r,t)p{r,t))e lEt / h dt 
-OO 

= 2irhN{r)n e (r,E), (102) 

wherein Af(r) is the mean density per unit volume of 
protons and n e (r,E) is the mean collective density per 
unit volume per unit energy of electrons which sit directly 
on the protons. 

The steady state inhomogeneous production of neu- 
trons per unit time per unit volume vj{t) as estimated 
in Eq. tjlOip ; i.e. exhibiting the radiated neutrino energy 



mc 2 f hG F 
irh \ c 3 



(gl- + 3<? A )AA(r)/C(r) 



JC(r) = - / n e {r. E = c 2 AM + huj)ujdLu, (103) 
c Jo 

wherein n e (r,E) must be calculated including the sur- 
face radiation energy and the driving current through 
the cathode. 

The calculation of K. depends on the detailed physical 
properties of the cathode surface as well as the flux $ 
per unit area per unit time of electrons determining the 
chemical cell current as in the power Eq.f3]). In the most 
simple model, let us consider a smooth surface with ma- 
terial properties and electron currents determining the 
neutron creation rate via the mass renormalization pa- 
rameter f3 as defined in Eq. (|48|) . We note in passing 
that a smooth surface is not likely to be the best sur- 
face for producing neutrons since rough surfaces have 
patches wherein the surface plasma electromagnetic field 
oscillations will be very much more intense. However, 
the following smooth surface model will be employed for 
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estimating the proper low energy nuclear reaction rates 
produced by electroweak interactions. 

To compute the density of surface electron states per 
unit area per unit energy, one may begin with a simple 
renormalized electron mass m model and take 

g 2 (E) = 2 j 5 (E - VcV + (mc 2 ) 2 ) , (104) 

i.e. 

92(E) = -§-2 . (105) 

If such surface electron states are confined to a wave func- 
tion width I normal to the surface, then we have within 
the surface region, and after integrating over the emitted 
neutrino energy spectrum 

/i f 

/C ~ — / gi(E — Pome 2 + hui)ujduj, 
Ic J 

^(f) 3 *"-*) 2 . (106) 

For a smooth surface, integrating the neutrino produc- 
tion rate over a thin slab at the electrode surface yields 
the estimate for the production rate per unit time per 
unit area, 

_ ( 9v + 39a \ ( G F m 2 \ 2 mc 2 2 

) n2 \—) (10?) 

The above Eq. (|107p for smooth surfaces is in agreement 
with the initial order of magnitude estimate in Eqs.@ 
and (HDJ. 

VII. CONCLUSIONS 

Electromagnetic surface plasma oscillation energies in 
hydrogen-loaded metal cathodes may be combined with 
the normal electron-proton rest mass energies to allow for 
neutron producing low energy nuclear reactions Eq.([2]). 
The entire process of neutron production near metallic 
hydride surfaces may be understood in terms of the stan- 
dard model for electroweak interactions. The produced 
neutrons have ultra low momentum since the wavelength 
is that of a low mode isotopic spin wave spanning a 
surface patch. The radiation energy required for such 
ultra low momentum neutron production may be sup- 
plied by the applied voltage required to push a strong 
charged current across the metallic hydride cathode sur- 
face. Alternatively, low energy nuclear reactions may be 



induced directly by laser radiation energy applied to a 
cathode surface. They may even be induced, albeit at 
comparatively low rates of reaction, simply by imposing 
an adequate pressure gradient (which can be as little as 
one atmosphere.) of gaseous hydrogen isotopes across 
a metallic membrane [32l l33j composed of an aggressive 
hydride-former such as palladium. 

The electroweak rates of the resulting ultra low mo- 
mentum neutron production are computed from the 
above considerations. In terms of the radiation induced 
mass renormalization parameter (3 in Eqs.® and 1(5)), the 
predicted neutron production rates per unit area per unit 
time have the form 

zu 2 = v(j3 - A)) 2 above threshold (3 > (3 . (108) 

The expected range of the parameter v for hydrogen- 
loaded cathodes is approximately 

10 12 < v < 10 14 (109) 
cm z cm 2 

in satisfactory agreement with the orders of magnitude 
observed experimentally. More precise theoretical es- 
timates of v require specific material science informa- 
tion about the physical state of cathode surfaces which 
must then be studied in detail. As discussed in previous 
workpj}, a deuteron on certain cathode surfaces may also 
capture an electron producing two ultra low momentum 
neutrons and a neutrino. The neutron production rates 
for heavy water systems are thereby somewhat enhanced. 

From a technological perspective, we note that energy 
must first be put into a given metallic hydride system in 
order to renormalize electron masses and reach the crit- 
ical threshold values at which neutron production can 
occur. Net excess energy, actually released and observed 
at the physical device level, is the result of a complex 
interplay between the percentage of total surface area 
having micron-scale E and B field strengths high enough 
to create neutrons and elemental isotopic composition of 
near-surface target nuclei exposed to local fluxes of read- 
ily captured ultra low momentum neutrons. In many 
respects, low temperature and pressure low energy nu- 
clear reactions in condensed matter systems resemble r- 
and s-process nucleosynthetic reactions in stars. Lastly, 
successful fabrication and operation of long lasting en- 
ergy producing devices with high percentages of nuclear 
active surface areas will require nanoscale control over 
surface composition, geometry and local field strengths. 
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